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ABSTRACT 



Most galaxy evolutionary models require quasar feedback to regulate star formation in their host galaxies. In particular, at high 
redshift, models expect that feedback associated with quasar-driven outflows is so efficient that the gas in the host galaxy is largely 
swept away or heated up, hence suppressing star formation in massive galaxies. We observationally investigate this phenomenon 
by using VLT-SINFONI integral field spectroscopy of the luminous quasar 2QZJ002830.4-281706 at z=2.4. The spectra sample the 
optical emission lines redshifted into the near-IR. The [OIII]/15007 emission-line kinematics map reveals a massive outflow on scales 
of several kpc. The detection of narrow Ho- emission reveals star formation in the quasar host galaxy, with SFR ~ 100 M Q yr~'. 
However, the star formation is not distributed uniformly, but is strongly suppressed in the region with the highest outflow velocity and 
highest velocity dispersion. This result indicates that star formation in this region is strongly quenched by the quasar outflow, which 
is cleaning the galaxy disk of its molecular gas. This is one of the first direct observational proofs of quasar feedback quenching the 
star formation at high redshift. 

Key words. Galaxies: formation - Galaxies: high-redshift - Galaxies: evolution - quasars: emission lines 



1. Introduction 

Most of the recent galaxy formation models invoke energetic 
outflows as a way to reg ulate the evolution of galaxies through- 
out the cosmic epoc hs dSilk & ReesL Il998t iBower et all I.2006; 
Springel et al., 2005]). In particular, quasars are expected to drive 
powerful outflows that eventually expel most of the gas in their 
host galaxies, thereby qu enching both star formation and fur- 
ther bl ack hole accretion dGranato et all l2004t p i Matteo et all 
l2005t IBower et all l2006t iHopk ins et al., 
e.g. ). According to those models, these 
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quasar driven outflows are required to prevent massive galaxies 
from overgrowing, hence explaining the shortage of very mas- 
sive galaxies in the local universe, and are responsible for the 
red color and gas poor properties of local elliptical galaxies. 

Massive, large-sc ale outflows have been detected in the hosts 
of local quasars (e . g. iFeruelio et all I2010I: iFischer et all [2010: 
ISturm et all 120111: iRupke & Veilleuxl 1201 ll) . HoweverTmodels 
expect that most of the quasar feedback action occurs at high red- 
shift, when quasars reach their peak activity (z~2) and when star 
formation in the most massive galaxies is observed to decline. 
Evidence of outflows in lu minous quasars has been detected 
up to very high redshift (e.g. lAllen et all 1201 it iMaiolino et all 
12004 lAlexander et all 1201 Oh . Indications that the strength of 
these AGN driven outflows anticorrelates with the starburst con- 
ttibution to the i nfrared luminosity has been obtained in red- 
dened quasars bv lFarrah et all d201 ll) . However, direct observa- 
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tional evidence that high-z quasar driven outflows quench star 
formation in their host galaxies is still missing. 

Here we present VLT-SINFONI near-IR integral field spectta 
of the quasar 2QZJ002830.4-281706 (hereafter 2QZ0028-28). 
This object, at z = 2.401, was take n from the sample of s trong 
[OIII]/15007 emitters discovered by Shem mer et all (|2004), and 
it is one of the most luminous quasars known. In this letter we 
show that the spatially resolved kinematics of the [OIII]/15007 
line clearly reveals a prominent outflow on scales of several 
kpc. Even more interesting, the star formation traced by nar- 
row Ha emission is suppressed in the region characterized by 
the strongest outflow. We suggest that this is the first observa- 
tional evidence of feedback associated with a quasar-driven out- 
flow that quenches star formation at high redshift, or one of the 
first. 



2. Observations and data reduction 

We used the near-IR integral field spectrometer SINFONI at 
the VLT to observe 2QZ0028-28 in the H and K bands, 
where H/3+[OIII]^4959,5007 and Ua+ [Nil] 6548, 65 84 are red- 
shifted, respectively, with the high-resolution gratings (deliv- 
ering R=3000 and R=4000, respectively). On-source exposure 
times were 1200s in the H band and 2400s in the K band (inter- 
leaved with sky observations). The observations were taken un- 
der excellent seeing conditions: 0.4" (measured from the broad 
lines, as discussed in the following). To properly sample this ex- 
cellent seeing, we used the camera that delivers a pixel scale of 
0".l x 0".05, providing a field of view of 3" x 3". 
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Fig. 2: Left panel: Velocity field (first moment map) of the [OIII]/15007 line, showing the prominent excess of blueshifted gas with 
a bow-like morphology SE of the nucleus. White contours are at 330, 360, and 390 km/s. Right panel: Velocity dispersion (second 
moment map) of the [OIII]/15007 line, showing the excess of dispersion in the SE region. White contours are at 680 and 700 km/s. 
In both maps the black contours trace the continuum. 




Fig. 1: Upper panel: 2QZ0028-28 H band spectrum extracted 
from the central 0.5 arcsec, along with the various components 
used for the fit (see Appendix [A] for details). Vertical dashed 
lines indicate the rest frame wavelength of each line, by taking 
the [OIII]/i5007 line peak for reference. Lower panel: Residuals 
of the fit. The green vertical lines enclose a the spectral zone 
affected by strong sky line residuals. 



Data was reduced with the ESO-SINFONI pipeline. The 
pipeline subtracts the background, performs the flat-fielding, 
spectrally calibrates each individual slice, and then reconstructs 
the cube. The pipeline delivers cubes where the spatial pixels are 
resampled to 0".05 x 0".05. 

3. Data analysis and results 

The final H and K band data cubes were analyzed by fitting each 
spatial spectrum separately in the field of view. Details of the 
fitting procedure for the spectra are given in AppendixlAl 



3.1. [OIII]+Hfi analysis: evidence for a powerful outflow 

The H-band spectrum (Fig. [TJ shows a clear broad H/3, associ- 
ated with the broad line region (BLR), and a prominent [OIII] 
doublet, primarily associated with the quasar narrow line region 
(NLR). The [OIII]/15007 line shows a clear asymmetric profile 
with a prominent blueshifted wing. This is generally a signature 
of outflows. Indeed, the lack of a corresponding redshifted wing 
is generally due to dust (in the host galaxy) absorbing the wind 
component on the opposite side of the galaxy disk, relative to 
our line of sight. 

The [OIII] line is well fitted with two Gaussians as seen in 
Fig. Q] However, since the two [OIII] components are heavily 
blended, their fits are not independent. To avoid the possible am- 
biguities associated with the two components possibly being cor- 
related, we investigate the [OIII] line kinematics by mapping the 
first and second moments of the global [OIII] profile (i.e. average 
velocity shift and dispersion). These two maps are shown in Fig. 
[2] Regions with S/N<3 were masked out. As a zero velocity ref- 
erence we have taken the peak of the [OIII] line in the central re- 
gion (Fig.[TJ. The black contours indicate the location of the con- 
tinuum. As shown by Fig.|2jleft), the average velocity is strongly 
blueshifted with respect to the line peak, by a few hundred km/s, 
as a consequence of the strong blue asymmetry of the line. The 
velocity is strongly negative over the whole region where [OIII] 
is detected, suggesting that, if the prominent blue wing is due 
to a wind, then the outflowing ionized gas is distributed over all 
of the central ~1 kpc. However, the most interesting features are 
the strongly blueshifted regions to the south and to the east of 
the nucleus, together making a bow-like morphology, suggestive 
of the envelope of a strong conical outflow. However, the ob- 
served velocity gradient may also be associated with rotation of 
the host galaxy disk (although the velocity field is certainly not 
a typical rotation curve). Outflow and galaxy rotation can be dis- 
tinguished through the velocity dispersion: outflows are always 
associated with high velocity dis persion compared with regu - 
lar motions in galaxy disks (e.g. Muller-Sanc hez et"aTl I201 ll) . 
The velocity dispersion map (Fig. [2j shows a large velocity dis- 
persion over the whole region where [OIII] is clearly detected, 
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suggesting that the whole central region is characterized by out- 
flows. However, the velocity dispersion increases significantly 
in the SE region, where the strongly blueshifted gas is detected. 
This correspondence between blueshift and velocity dispersion 
strongly suggests that the blueshift observed in this region is due 
to an outflow excess and not to galaxy rotation. 





H %J 










Ha [Nil] 

:.ltiJ..yjri, Jj ..„ J lJ.„ 













Wavelength [f+m] 

Fig. 3: Top panel: 2QZ0028-28 K-band spectrum extracted from 
the central arcsec, along with the various components used for 
the fit (blue are for Ha components, magenta are for [Nil] lines 
of component B). Vertical dotted lines show the rest-frame wave- 
length for Ha and [NII]/t6584, by using the same reference as in 
Fig. [TJ Middle panel: Residuals of the fit. Bottom panel: Result 
of the subtraction between the spectrum extracted from the re- 
gion with narrow Ha emission NW of the nucleus (Fig. |4]i and 
the spectrum extracted from the region without narrow Ha emis- 
sion, SE of the nucleus, after scaling the two spectra to match the 
intensity of the broad line. A clear narrow Ha component is de- 
tected, illustrating that the detection of this component and its 
distribution are not artifacts of the spectral fit. 



driven outflows dThackeret all [2006). The mass of ionized out- 
flowing gas is simply given by (see AppendixlBl) 



MT = 5 3 10 7 L44([ ° m]) 



(1) 



where L 44 ([OIII]) is the luminosity of the [OIII]/15007 line trac- 
ing the outflow, in units of 10 44 erg/s, n e3 is the electron density 
in the outflowing gas, in units of 10 3 cm 3 , typical of the NLR, 
and 10 [O/,H1 is the oxygen abundance is Solar units. If we assume 
that the outflow is traced by the broad component of [OIII], then 
L44([OIII]) ~ 2, implying that the mass of ionized gas involved 
in the outflow is ~ 10 8 M . 

By assuming a simplified conical (or biconical) outflow dis- 
tributed out to a radius Rk pc (in units of kpc), then the outflow 
rate of ionized gas is given by (see AppendixlBl 



M? 



164 



L 4 4([OIII]) v 3 
n e3 IOIO/h] R kpc 



M yr 



(2) 



where v 3 is the outflow velocity in units of 1000 km s . 
The maximum outflow velocity inferred from the [OIII] profile 
is about 2000 km/s, which is probably representative of the av- 
erage outflow velocity, while the lower velocities observed in 
the [OIII] line profile are very likely projection effects. By us- 
ing Rkp C = 3, we infer an ionized gas outflow rate of about 
200 M Q yr . This is a lower limit of the total outflow rate, since 
the ionized component is probably a minor fraction of the global 
outflowing gas mass. If we scale by the same neutral-to-ionized 
fraction as in Mrk231 (Rupke, priv. comm.), the total outflow 
rate can be up to an order of magnitude higher. 

The inferred kinetic power of the ionized component of the 
outflow is about 3 10 44 erg s _1 (Eq ]B.9l >, or about ~ 10~ 3 times 
the bolometric luminosity of the AGNs (~ 3.8 10 46 erg s~ l , from 
the continuum luminosity at 5100A, and by assuming a bolo- 
metric correction factor of 7). However, if the ionized outflow is 
accompaneid by a neutral/molecular outflow an order of magni- 
tude more massive, as discussed above, then the kinetic power 
is also likely to be an order of magnitude higher, i.e. about ~1% 
of the bolometric luminosity , close to the expectation of quasar 
feedback models dLapi et al.ll2005l) . 



While there is a general correspondence between velocity 
blueshift and dispersion, both with the highest values in the SE 
quadrant, a detailed analysis reveals interesting differences. The 
peak of the velocity dispersion is located between the two peaks 
of velocity blueshift. A possible interpretation is that those ar- 
eas with highest velocity dispersion are associated with regions 
where the outflow is strongly interacting with the gas in the host 
galaxy disk, which also lowers the velocity locally. 

The outflow may even be intrinsically symmetric, but 
the opposite (receding) outflow is likely to be obscured by 
dust in the host galaxy d isk (as observed in local AGNs, 
Mull er-Sanchez et all 1201 II) . It is interesting to note that the ve- 
locity dispersion map shows a strong excess in a tiny region on 
the NW (opposite to the main [OIII] outflow), this may be trac- 
ing the opposite outflow coming out of the galactic dusty disk. 

The [OIII] line is primarily tracing gas in the NLR ionized 
by the QSO. Therefore, the most likely explanation is that the 
outflow is driven by the QSO radiation pressure. A quasar ori- 
gin of the wind is also supported by the fact that the outflow- 
ing gas reaches velocities in excess of 1000 km/s (in the [OIII] 
blue wing), which cannot be explained by models of supernovae- 



3.2. Ha+[NII] analysis: evidence for quenched star formation 

The Ha has a very broad profile (from the BLR), and its over- 
all profile has been fitted mostly with two broad Gaussians (see 
Appendix lAl for details). Interestingly, the fit also requires the 
presence of a weak, but significant narrow component of Ha 
(A*, FWHM~ 600/bra/s). The absence of [Nil] indicates that 
this narrow Ha emission is mostly due to star formation. 

The map of the narrow Ha emission is shown in Figj4] which 
reveals star formation extending over a few kpc from the nu- 
cleus. However, the star formation is not distributed symmetri- 
cally around the nucleus, but primarily towards the N and W. In 
particular, the region within a few kpc SE of the nucleus is nearly 
free of any narrow Ha tracing star formation. To show that the 
detection of narrow Ha and that the asymmetric distribution are 
not artifacts of the fitting procedure, we extracted a spectrum by 
integrating over the NW region putatively containing star for- 
mation, according to the Ha map, and another spectrum in the 
SE region devoid of star formation. After scaling the two spec- 
tra so that the intensity of the broad Ha wings (associated with 
the BLR) are the same, we have obtained the difference of the 
two spectra, which is shown in the bottom panel of Fig. [3] The 
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Fig. 4: Map of the narrow component of Ha with contours tracing the [OIII] velocity shift (left panel) and velocity dispersion (right 
panel), as in Fig. [2] Star formation, traced by Ha, is heavily suppressed in the SE region where the strongest outflow is traced by 
[OIII]. 



differential spectrum clearly shows the presence of the narrow 
Ha. The line broadened profile (also making the detection nois- 
ier than expected ) is because we are integrating over different 
regions of the velocity field. The differential spectrum also con- 
firms that no [Nil] is detected in association with the narrow 
Ha, confirming that the latter is tracing star formation and not 
the quasar NLR. 

We mention that the map of the [OIII] narrow component 
(shown in the appendix) is also characterized by a similar asym- 
metry towards the NE, suggesting that some of the [OIII] narrow 
line is associated with star formation. However, the asymmetry 
is less clean than observed for the Ha narrow component, likely 
because of NLR contribution to [OIII] and because of the blend- 
ing with the "broad" [OIII] component. 

The integrated emission of the narrow Ha yields a total star 
formation rate in the host galax y of about 100 M B yr" 1 (by using 
the conversion factor given in Kennicutt, 1998), which is not 
unusual in high-z quasars (e.g. iLutz et all 12008!). However, the 
most interesting result is that the star formation is heavily sup- 
pressed in the SE region, which is characterized by the excess of 
outflow with high-velocity dispersion. In Fig. |4|left) the white 
contours identify the strongest gas outflow traced by the highly 
blueshifted [OIII] line, as in Fig.|2]-left, while in Fig.HJ-right the 
white contours identify the highest velocity dispersion region, 
as in Fig. fright, which is likely the region where the strong 
outflow interacts with the host galaxy disk. We suggest that the 
heavy suppression of star formation in the region of strongest 
quasar-driven outflow among the first direct observational proofs 
of quasar feedback onto the host galaxy quenching star forma- 
tion at high redshift, as predicted by models. 

4. Conclusions 

By using near-IR integral field spectroscopic observations we 
have revealed a powerful outflow in the host galaxy of the quasar 
2QZ0028-28 at z=2.4. The outflow was revealed by the veloc- 
ity field traced by the [OIII]/i5007 line, redshifted into the H- 
band. We estimated that the outflow rate of ionized gas is about 
200 M yr~', which is, however, a lower limit of the total gas 



outflow rate. Both the high outflow velocity (> 1000 km/s) and 
the fact that the wind is mostly traced by the [OIII] line (pro- 
duced primarily in the NLR) strongly suggest that the outflow is 
mostly driven by the quasar. The outflow is not symmetric, the 
highest velocities and highest velocity dispersion are found in 
the region SE of the nucleus. 

In the K-band, our data clearly reveal the presence of narrow 
Ha emission tracing star formation in the host galaxy, on scales 
of several kpc and with a rate of about 100 M yr _1 . However, 
star formation is not distributed uniformly in the host galaxy, 
but is mostly found in the regions not directly invested by the 
strong outflow. Instead, star formation is heavily suppressed in 
the SE region where the strongest outflow is detected. This ob- 
servational result supports models invoking quasar feedback to 
quench star formation in massive galaxies at high redshift. 
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Appendix A: Details of the spectral fitting 

The initial spectral fit is performed on a spectrum extracted from 
a central aperture of ten pixels (i.e. 0.5 arcsec) in the H band and 
20 pixels (i.e. 1.0 arcsec) in the K band, which guarantee high 
S/N. 

The emission lines are fitted with multiple Gaussians and, 
in th e case of the broad lines, by also using power-law profiles, 
as in Nagao et al. (2006). The continuum is fitted with a single 
power law, which represents the emission of the QSO accretion 
disk (plus possibly some minor contribution from the host galaxy 
stellar continuum). Starting from this initial fit in the central re- 
gion, we fitted the spectra individually at all spatial pixels by 
leaving most of the parameters free, except for the velocity, dis- 
persion, and relative intensity of the components describing the 
broad lines (Ha and H/3). Indeed, since the BLR is unresolved, 
the shape and shift of the broad line profile (components C+D+E 
in the case of H/3) must be constant over the field of view, and 
there the global intensity variation only reflects the seeing PSF. 
We tri ed to also inc l ude an Fe II template in the spectral fitting 
(as in iNetzer et all 120041) : however, this is always set to zero 
by the fitting procedure, confirming the lack of significant Fell 
emission inferred by the visual inspection of the spectrum. Thi s 
result contrasts with the measurement of INetzer et al.l (|2004), 
who obtain a flux of the Fe II emission that is about 0.37 times 
the H/3 emission, and we ascribe the discrepancy to the much 
lower signal-to-noise in the latter spectrum. 

Figure [TJ shows the various components used to fit the H- 
band spectrum of the central region and, in the bottom panel, 
the fit residuals. The region included within the two green lines 
is affected by strong sky residuals and was not considered in 
the fit. In this figure we can see that the [OIII] profile can 
be nicely fitted with two Gaussians, one relatively narrow (A, 
FWHM~600 km/s) and a second one (B) blueshifted by about 
700 km/s and with FWHM-1700 km/s (much broader than 
typically found in the NLR of lower luminosity AGNs). The 
[OIII]4959 line is fitted with the same components, linked to 
have an intensity equal to one third of the [ODTJ/15007 line. 
Figure I A. 1 1 shows the flux distribution of the broad and narrow 
components of [OIII]. The peak of the [OIII] broad component 
is shifted towards the SE, i.e. in the same direction as the strong 
outflow, confirming that this is the region where the NLR de- 
velops and where the quasar radiation pressure is driving the 
outflow. The narrow component of [OIII] is instead distributed 
towards the N and towards the W, i.e. similar to the Ha narrow 
component, suggesting that the narrow component of [OIII] re- 
ceives a significant contribution from star formation. However, 
the imperfect correspondence between the two maps suggests 
that a fraction of the narrow [OIII] is also contributed by the 
NLR. 

The [OIII]4959 line is heavily blended with another line at 
^rest - 4930A, which is also seen in the spectra of other quasars, 
sometimes tentatively identified with Fell emission; however, in 



the quasar discussed here (as well as in a few others showing 
the same feature), the iron emission that typically encompasses 
the H/3+[OIII] group is particularly weak, therefore suggesting a 
different origin of this line. In the residuals we observe a narrow 
component that is likely associated to the same unidentified line. 
The fit of the broad Hf3 requires three Gaussians and a powerlaw 
profile (as in Nagao et al. 2006). We also included the H/3 contri- 
bution associated with the two [OIII] components. However, one 
should bear in mind that the intensity of these weaker H/3 com- 
ponents is difficult to evaluate, since these are blended within 
the broad, complex H/3 profile. In Table lA.fl we list the parame- 
ters inferred for all of the components used in the fitting of the 
SINFONI spectra. 

The Ho- profile (dominated by the broad component, trac- 
ing the BLR) is clearly different from the H/3 profile, which is a 
property that is common to many other quasars and AGNs and 
that is ascribed to complex radiative transfer within the dense 
gas of the BLR. The bulk of the Ha profile was fitted by using 
two very broad Gaussians (H and G in Fig. |3), which give the 
seeing in the K-band, which is roughly consistent with the see- 
ing measured in the H-band observation. As for the broad H/3, 
the relative intensity, shift, and width of these two lines are kept 
fixed over the field of view, and their overall intensity variation 
reflects the seeing PSF. 

We forced the inclusion of component B in the Ha profile, 
by using the same velocity shift and width as the correspond- 
ing [OIII] component, while the intensity was left free to vary. 
This component of Ha is certainly associated with the NLR. 
However, since this component is relatively broad, its intensity 
is poorly constrained (see uncertainty in the flux of this com- 
ponent in Table IA.ll . and degenerate with the other two broad 
Ha Gaussian components (as is the case for the corresponding 
component in H/3). The [Nil] doublet associated with compo- 
nent B is also forced to have the same shift and width as the 
corresponding [OIII] component. The relative intensity of the 
two [NII]6548,6584 lines are forced to be in the ratio of 1:3. 
The intensity of these [Nil] lines is even less constrained than 
the corresponding Ha B component, since the wavelength of 
[NII]6548 nearly overlaps the intense component G of Ha. The 
resulting [NII]/Ha ratio of component B is low, but highly uncer- 
tain (log (F[pfji]/Ffi a ) = -0.73 + 0.45), but still consistent with 
the values observed in AGNs. 

The narrow component "A" is much narrower than the other 
components and, therefore, easier to disentangle from the broad 
Ha profile. However, in principle, we should include two nar- 
row (FWHM ~ 600km/s) Ha components, one associated with 
the quasar NLR and another one associated with any putative 
star formation in the host galaxy. However, the quality of our 
data does not really allow us to fit two separate narrow compo- 
nents, because these would be totally degenerate. We therefore 
fit a single narrow component not tied to have the same velocity 
and width of component "A" of [OIII]. We label this narrow Ha 
component with "A*" (meaning that it may be partly associated 
with component A of [OIII], but not necessarily). We investigate 
the relation of the Ha component "A*" with the [OIII] compo- 
nent "A" a posteriori. We note that it is not possible to follow 
a similar approach on H/3 (i.e. introduce a component "A*", not 
linked to the [OIII] components, since, besides the problem of 
the blending with other components, the signal-to-noise on H/3 
is much lower). 

The parameters resulting from the best fit in the K-band are 
given in Table lA.ll We note that there is no room for narrow A* 
[Nil] emission, at a level of F[nh]6584(A*) < 0.21 FH ff (A*). The 
lack of an [Nil] narrow component is confirmed by the "differ- 
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Fig. A.l: Flux maps of the broad (B, left) and narrow (A, right) components of [OIII]. 



ential" spectrum presented in Sect. 13.21 and in Fig. [3] which is 
totally independent of any fitting procedure. As discussed in the 
body of the paper, the lack of [Nil] at a level below one fourth 
of Ho; indicates that this narrow Ha emission is mostly tracing 
star formation, and not the NLR. 



Figure IAT2I shows the velocity field and the velocity disper- 
sion of the narrow component of Ha. Both maps are very noisy, 
owing to the weakness of the line. The velocity field does not 
clearly indicate a rotation pattern, which would be expected by 
gas in a galactic disk, except possibly for an NW-SE gradient, 
but the latter may be associated with some contribution to Ha 
narrow from the outflow in the SE region. However, only a frac- 
tion of the disk is actually traced by the Ha narrow, and this, 
together with the noisy velocity map, may prevent identification 
of a clear rotation pattern. Moreover, it is well known that host 
galaxy disks of optically selected quasars tend to be face on, as 
a consequence of selection effects (Carilli & WangL 12006). so it 
is not expected that quasar host galaxies have prominent rota- 
tional patterns. The velocity dispersion map is very noisy, but it 
is consistent with being uniform over the area where Ha narrow 
is detected. 

We finally note that the best-fit velocity and FWHM of com- 
ponent A* of Ha are similar to component A of [OIII]. This 
further suggests that the latter component of [OIII] is partly con- 
tributed by the ionized gas in the star-forming regions traced by 
the narrow Ha. The inferred F^'/F^m ~ 0.17, at the verge of 
the range typically observed in star-forming galaxies, indicates 
that the flux of component A of [OIII] is not incompatible with 
being partly originated by star formation, but probably a contri- 
bution by the AGN NLR is required. As mentioned above, the 
similarity of the F[omj(A) map and the Fhq.(A*) map also sup- 
ports the scenario where part of component A of [OIII] is asso- 
ciated with star formation. 



Appendix B: A simple model of the ionized outflow 

In this section we discuss how the physical properties of the ion- 
ized outflow can be constrained through the observational pa- 
rameters of the [OIII] line, by adopting a simple model for the 



ionized wind. The [OIII]5007 line luminosity associated with the 
outflow is simply given by 



L([OIII]) 



qom] f dV 



(B.l) 



where V is the volume occupied by the outflowing ionized gas, f 
the filling factor of the [OIII] emitting clouds in the outflow, and 
f[Oiii] the [OIII]5007 emissivity that, at the temperature typical of 
the NLR (~ 10 4 K), has a weak dependence on the temperature 
(oc T 01 ). It can be expressed by 



qom] = 1-11 10 hv [0 iii] n 0+ 2 n e erg s cm 



(B.2) 



where hv\pm\ is the energy of the [OIII]5007 photons (in units of 
erg), no+2 and n e are the volume densities of the +2 ions and of 
electrons, respectively (in units of cm -3 ). Under the reasonable 
assumption that most of the oxygen in the ionized outflow is in 
the +2 form, then 



qom] * 5 10- 13 hv [0 irn n e 2 10 l ° /H] erg s -1 crrT 3 

where 10 [O ^ H1 gives the oxygen abundance in solar units. 
The mass of outflowing ionized gas is given by 



MH 



X 



1.27m H n e f dV 



(B.3) 



(B.4) 



where mn is the mass of the hydrogen atom, and where we have 
neglected the mass contributed by species heavier than helium. 
By combining Eqs. IB.ll and lB.4l we obtain 



(B.5) 



M io„_ 5 33 10 7 CL44([oni]) 

M out - 5.33 10 <ne3> 1Q[0/H] M 



where L 44 ([OIII]) is the luminosity of the [OIII]5007 line 
emitted by the outflow, in units of 10 44 erg s _1 , (n e 3) (= 
J v n e f dV/ J v f dV) is the average electron density in the ion- 
ized gas clouds, in units of 10 3 cirT 3 , and C = (n e 3) 2 /(n 2 3 ) is a 
"condensation factor", where (n 2 3 ) = J n 2 f dV/ J f dV. We can 
assume C = 1 under the simplifying hypothesis that all ionizing 
gas clouds have the same density. Also, under these assumptions, 
the mass of outflowing ionized gas is independent of the filling 
factor of the emitting clouds. 
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Fig. A. 2: Velocity field (left) and velocity dispersion (right) of the narrow component (A*) of Hot. 



Line 


Component 


Fitting Function 




FWHM 


Flux 


Velocity 1 








w 


(km/s) 


(lO^ergcnT 2 s" 1 ) 


(km/s) 


[OIII]5007 


A 


Gaussian 


1.7061 


652 


10.9+0.56 


92 




B 


Gaussian 


1.7020 


1797 


43.6+2.3 


-663 


w 


A 


Gaussian 


1.6564 


652 


0.79+0.66 


92 




B 


Gaussian 


1.6522 


1797 


8.8±3.6 


-663 




C 


Broken Power Law 


1.6567 


5508 


57.9±4.6 


156 




D 


Gaussian 


1.6648 


1150 


4.4±2.9 


1622 




E 


Gaussian 


1.6431 


2243 


11.3+0.26 


-2306.6 


No ID 


F 


Gaussian 


1.6806 


2378 


13.18+2.3 




Ha 


A* 


Gaussian 


2.2361 


616 


1.9+0.3 


20 




B 


Gaussian 


2.2310 


1797 


9.3±3.5 


-663 




G 


Gaussian 


2.2391 


3386 


73.7+5.5 


422 




H 


Gaussian 


2.2478 


9717 


220.8±10.4 


1588 


[NII]6584 


A* 


Gaussian 


2.2433 1 


616' 


<0.4 


I 1 




B 


Gaussian 


2.2382 


1797 


1.7+1.5 


-663 



Table A.l: Results of the spectral fit for the individual components. 

The names of the components correspond to those shown in Figs. Q] and [3] Notes: 1 The velocities of the components are obtained 
by assuming the peak of [OIII]/15007 as a rest frame reference. 2 For this [Nil] line, the width and velocity were forced, to estimate 
the upper limit to the values found for the corresponding A* component in Ha. 



If we assume a simplified model of the outflow (justified by 
the limited information currently available to us) where the wind 
occurs in a conical region, with opening angle Q., composed of 
ionized clouds uniformly distributed and outflowing with veloc- 
ity v, out to a radius R, then the mass outflow rate of ionized gas 
is given by 



MI! 



<p)v v QR 



(B.6) 



where (p)y is the average mass density in the whole volume oc- 
cupied by the outflow, which is given by 



jyrion 



(B.7) 



where the volume occupied by the conical outflow is given by 



V 



±7TR 



3 £2 



,„n i . Unless f = 1, generally (p)y + 1.27 mn (n e ), since 
the latter numerical density (defined above) is averaged among 
the emitting clouds, not over the whole volume. 



By replacing Eqs. IB. 5 1 and IB. 71 into Eq. IB. 61 we obtain that 
the ionized outflow rate is given by 



MI! 



164 



CL44([OIH])v 3 
<n e3 ) 10[°/H] R kpc 



M yr 



(B.8) 



where L44([OIII]), n e 3, and C (~ 1) were defined above, V3 is 
the outflow velocity in units of 1000 km s _1 , and Rk pc is the 
radius of the outflowing region, in units of kpc. The outflow rate 
is independent of both the opening angle Q of the outflow and of 
the filling factor f of the emitting clouds (under the assumption 
of clouds with the same density). 

The kinetic power (associated with the ionized component) 
is then given by 



P£ n = 5.17 10 



,43 



CL 44 ([OIII])v3 
<n e3 ) lOlo/Hl R kpc 



ergs 



(B.9) 
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